Boundary-induced Anderson localization in photonic lattices 
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We analyze numerically localization of light in linear square waveguide arrays restricted in one 
dimension ("ribbons"), whose boundaries are disordered in propagation constant and/or coupling. 
We find that the disordered boundary induces a localization tendency in the bulk even for relatively 
wide ribbons. 
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The phenomenon of Anderson localization (AL) in dis- 
ordered system constitutes one of the staples of modern 
condensed-matter physics. Proposed originally for elec- 
trons and one-particle excitations in solids QUI], it was 
soon extended to many other fields such as acustics[5- 
[7], Bose-Einstein condensates [8] and optics [9HT3]. It still 
continues to surprise us even now, more than 50 years 
since its discovery [14]. 

The AL effect is ultimately due to wave interference be- 
tween multiple-scattering paths, which makes optical sys- 
tems an ideal setting for exploring the AL phenomenon, 
without encumbering effects such as electron-electron in- 
teractions in solids, or other many-body effects. Thus 
far, most research carried on in optical disordered sys- 
tems has focussed mainly on bulk disorder. However, in- 
teresting localization effects due to surface disorder have 
been found in single-and-multimode waveguides. For in- 
stance, localization of waves due to multiple-scattering by 
surface roughness has been predicted for a single-mode 
waveguide and for a thin film[15l-IT7]. 

In the same spirit, we consider in this work the case of a 
two-dimensional photonic lattice with boundary disorder. 
As we show below, for finite width ribbons, boundary 
disorder does induce AL in the bulk. 

Let us consider a two-dimensional square linear N x M 
waveguide array with a finite extension in one dimension 
(Fig.l), i.e., a "ribbon" (M <C N). In the framework 
of the coupled-modes theory, the electric field E(r, z) 
propagating along the waveguides can be presented as 
a superposition of the waveguide modes, E(r, z) = 
J2 n En(z)<fi(r — n), where r = (x,y), E n is the ampli- 
tude (in units of (W) 1 / 2 ) of the single guide mode 0(r) 
centered on site with the lattice number n = (774,77,2). 
The evolution equations for the modal amplitudes E n 
are 

i^^ + e n E n (z) + Y,Vn, m E m (z)=0, (1) 

m 

where n denotes the position of the guide center, z is the 
longitudinal distance (in meters), T4i,m is the coupling 
between nearest-neighbor guides (in units of 1/m) and 
e n is the propagation constant of guide with center at n. 

Stationary states. The stationary modes of the sys- 
tem are found by posing a solution of the form E n (z) = 
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FIG. 1: Example of optical lattice ribbon (M <C N) with 
disordered boundary. The propagation constants along the 
long boundary have random values. 

C n exp(if3z) in EqQ, which leads to 

- (3C n + e n C n + ^n, m C m = 0. (2) 

m 

To ascertain the localization properties of a given mode of 
the system, we use the inverse participation ratio (IPR), 
defined by IPR = £ n |C n | 4 /^ n |C n | 2 . For completely 
localized modes, IPR = 1, while for completely delocal- 
ized modes, IPR = 1/N, where N is the number of sites 
of the lattice. We will average the IPR over all states 
of the system and also over a number of random realiza- 
tions, to obtain the average IPR which reflects the global 
localization tendency of the system. We start with the 
case where the propagation constants e n along the bound- 
ary of the ribbon (gray sites on Fig{T]) take on random 
values stemming from a uniform distribution of width w: 
e n G [—w : w] : while coupling between nearest-neighbor 
guides is identical (= 1). Results for the average IPR are 
shown on Figj2|a), for several disorder widths, including 
the case of no disorder w = 0, which is necessary for com- 
parison, since we are dealing with finite systems. Figure 
|2|b) shows the corresponding results for a ribbon with all 
its propagation constants random. As expected, the IPR 
is smaller for the surface disorder case than for the bulk 
disorder case. But we can also see that for boundary dis- 
order, some localization tendency is still appreciable for 
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FIG. 2: State-and-realization average inverse participation ra- 
tio vs ribbon width for (a)-(b) disorder in waveguide propaga- 
tion constant, and (c)-(d) disorder in waveguide coupling, for 
several disorder widths: w = O(circles), w = 1 (triangles) and 
w = 2 (squares). Left (right) column corresponds to surface 
(bulk) disorder. 



intermediate widths, e.g., M = 4, 5. For the case where 
e n = 0, but the couplings along the boundary of the 
ribbon (between shaded sites on Fig{T]) take on random 
values from a uniform distribution: V n , m = 1+5V, where 
SV G [—w/4,w/4], Fig{2] shows qualitatively similar re- 
sults as for the disordered propagation constants case, al- 
though the localization effect is smaller in this case. The 
above results suggest the idea of a disorder "memory" in 
intermediate width boundary-disordered ribbons. 

Dynamics. We examine now the evolution (Eq{T]) of 
an initially localized input beam launched at the mid- 
dle site belonging to one of the not-disordered bound- 
aries of the ribbon, (0, M/2). The long boundaries are 
disordered, with random propagation constants e n . The 
idea is to examine the evolution of the spatial optical 
power distribution, and compare it to the cases with 
no disorder (e n = constant) and to the case with bulk 
disorder (e n = random for all sites). Since the prob- 
lem is linear, we can express E n (z) as a superposition of 
the stationary modes of the system. For localized initial 
conditions E n (0) = £ nno , with n = (0, M/2), the solu- 
tion of Eq.Q is E n (z) = Cj^C^e-W™ where c£° 
and /3(0 are the 1-th mode and its associated propagation 
constant of the eigenvalue Eq. Q . For the completely or- 
dered array (Fig{3|, the beam spreads out, away from the 
boundary, and most of its power is carried in the frontal 
lobes. In the case of bulk disorder (FigQ we see the 
usual Anderson localization phenomena with the optical 




FIG. 3: Dynamical evolution of an initially localized in- 
put beam at the boundary of a 41x5 completely ordered 
photonic ribbon. From left to right we show snapshots at 
Vz = 0,2,4,6,8,10,12,14,16 



power remaining confined within a small region around 
the excitation point. For the case of boundary disorder 




FIG. 4: Average over 100 realizations of dynamical evolu- 
tion of an initially localized input beam at the boundary 
of a 41x5 bulk-disordered (e n G [—1,1] for all n) photonic 
ribbon. From left to right we show snapshots at Vz = 
0,1,2,4,8,16,32,64,128. 

, however, something very interesting happens: After a 
short transient longitudinal distance, the extent of which 
depends on both, the width of the disorder w and the 
relative width of the ribbon M, the system tends to lo- 
calize a finite fraction of its power at the initial launching 
site, while the rest of the optical power gets eventually 
Anderson-localized (Figs {5] and [6|. The effect is stronger 
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FIG. 5: Average over 100 realizations of the dynamical evo- 
lution of an initially localized input beam at the boundary 
of a 41x5 boundary-disordered (e n G [—1,1] along the long 
boundaries) photonic ribbon. From left to right we show snap- 
shots at Vz = 0, 1, 2, 4, 8, 16, 32, 64, 128 

for narrow ribbons (e.g., 5 sites), but has been verified 
for wider ribbons as well, albeit, with a greater localiza- 
tion evolution distance ("time"). This is only natural, 
since for wider ribbons it takes longer to achieve enough 
multiple scattering for localization effects. For the case 
when the disordered lies only at the narrow edge where 
the beam is launched (narrow lower edge in Fig.l), the 
above localization effect ceases, and we only obtain par- 
tial localization at the input site, plus free propagation 
of the remaining fraction (not shown). Thus, it seems 
that multiple wave scattering between at least two disor- 
dered boundaries is needed to effect AL in these kind of 



photonic lattices. 
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FIG. 6: Time and realization (100) average power remaining 
at initial guide vs ribbon width, for disorder width w = 
(rhombi), w = 0.5 (squares), w = 1 (triangles) and w = 2 
(circles). 

Therefore, the presence of boundary disorder does ef- 
fect Anderson localization on the whole system. This ef- 
fect is reminiscent of localization in a single- mode waveg- 
uide with rough edges [15 , and suppression of thermal 
conductivity in thin graphene nanoribbons with rough 
edges [18] . 

In summary, we have analyzed numerically localiza- 
tion effects on a two-dimensional photonic lattice of finite 
width with disordered boundaries, by computing the av- 
erage inverse participation ratio of the stationary modes, 
and by dynamical evolution of an initially localized in- 
put beam. We conclude that a disordered boundary does 
induce AL effects on the whole bulk. 
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